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Lithium  ion  battery  electrode  materials  generally  experience  significant  volume  changes  during  charging 
and  discharging  caused  by  concentration  changes  within  the  host  particles.  Electrode  failure,  in  the  form 
of  fracture  or  decrepitation,  may  occur  as  a  result  of  a  highly  localized  stress,  strain  energy,  and  stress 
cycles  over  time.  In  this  paper,  we  develop  analytic  expressions  for  the  evolution  of  stress  and  strain  energy 
within  a  spherically  shaped  electrode  element  under  either  galvanostatic  (constant  current)  or  poten¬ 
tiostatic  (constant  potential)  operation  when  irreversible  phenomena  are  dominated  by  solute  diffusion 
resistance  within  host  particles.  We  show  that  stresses  and  strain  energy  can  evolve  quite  differently 
under  potentiostatic  vs.  galvanostatic  control.  The  findings  of  this  work  suggest  the  possibility  of  devel¬ 
oping  new  battery  charging  strategies  that  minimize  stress  and  strain  energy  and  thus  prolong  battery 
life. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Diffusion-induced  stress  (DIS)  can  occur  as  a  result  of  com¬ 
positional  inhomogeneities  during  solid-state  diffusion  in  many 
technological  situations,  including  dopant  diffusion  in  semicon¬ 
ductor  processing,  oxidation  of  metals,  hydrogen  transport  in 
solid-state  hydrogen-storage  media,  and  lithium  diffusion  in  bat¬ 
tery  electrodes.  Several  authors  have  developed  models  for  DIS. 
For  example,  Prussin  [1]  made  an  analogy  between  thermal  stress 
and  DIS  and  analyzed  the  transverse  stresses  developed  in  a  thin 
plate  during  mass  transfer.  Li  [2]  provided  a  number  of  analyt¬ 
ical  solutions  to  DIS  problems  in  spherical,  cylindrical,  and  thin 
plate  geometry.  Lee  and  co-workers  [3-5]  have  also  studied  DIS  in 
various  systems  including  thin  plates,  hollow  cylinders,  and  com¬ 
posites.  Yang  and  Li  [6,7]  considered  the  effect  of  diffusion-induced 
stresses  on  the  bending  of  beam  and  plate  structures  for  sensing 
applications.  Yang  also  studied  the  coupled  problem  of  interactions 
of  stress  and  diffusion  [7]. 

Within  the  electrochemical  energy  storage  topical  area,  a  num¬ 
ber  of  publications  address  the  modeling  of  lithium  ion  batteries 
[8-15]  and  are  focused  on  the  underlying  thermodynamics,  interfa¬ 
cial  kinetics,  and  transport  phenomena  governing  electrochemical 
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systems  in  the  context  of  a  volume-averaged  treatment  of  porous 
media.  Recently,  Garcia  et  al.  [16],  Christensen  and  Newman  [17,18] 
and  Zhang  et  al.  [19]  studied  stresses  generated  during  Li  diffusion 
under  galvanostatic  control.  The  latter  authors  also  investigated 
linear  sweep  voltammetry  of  a  single  particle  as  well  as  heat  gen¬ 
eration  [20];  their  investigation  included  spherical  and  ellipsoidal 
particles.  The  paper  by  Zhang  et  al.  [21]  is  devoted  to  the  analy¬ 
sis  of  a  single  particle  undergoing  electrochemical  insertion  and 
outlines  the  motivation  for  such  single  particle  studies.  The  afore¬ 
mentioned  battery  references  and  this  work  are  motivated  by  the 
need  to  develop  high  energy  and  power  batteries  that  are  durable 
over  the  intended  usage  profiles  [24]. 

Since  batteries  are  usually  charged  using  sophisticated  meth¬ 
ods  that  consist  of  both  galvanostatic  and  potentiostatic  control, 
we  report  in  this  paper  stress  and  strain  energy  evolution  under 
either  purely  galvanostatic  or  potentiostatic  condition  when  the 
resistance  of  the  cell  is  governed  by  solute  diffusion  limitations 
within  the  insertion  particles.  The  analytic  solutions  developed  in 
this  work  should  improve  the  current  understanding  of  stress  and 
strain  energy  evaluation  by  (1 )  identifying  dimensionless  parame¬ 
ters  that  control  stress  and  strain  energy  evolution,  (2)  providing 
the  asymptotic  behavior  for  short-  and  long-time  limits,  and  (3) 
rendering  order  of  magnitude  estimates  for  more  complicated  oper¬ 
ating  conditions.  The  results  of  this  work  may  be  used  to  develop 
strategies  that  minimize  stress-induced  battery  failures.  Additional 
information  on  relevant  past  works  by  the  present  authors  can  be 
found  in  Refs.  [22,23];  in  these  two  investigations,  surface  tension 
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Nomenclature 

List  of  symbols 

as 

specific  surface  (e.g.,  cm2  cm-3) 

C 

molar  concentration 

Co 

initial  (uniform)  solute  concentration 

Cr 

surface  concentration  during  constant  potential 
operation 

Qzv(0 

average  concentration  in  the  spherical  volume  of 
radius  r 

Cav(R) 

average  concentration  in  the  spherical  particle  of 
radius  R 

D 

diffusion  coefficient  of  the  solute 

E 

Young’s  modulus 

Et 

total  elastic  energy  stored  in  the  sphere  of  radius  R 

e 

strain  energy  density 

e~ 

electron 

E 

Faraday’s  constant 

H+ 

proton 

/ 

current  over  a  spherical  particle  of  radius  R 

hell 

cell  current 

h,ref 

exchange  current  based  on  reference  concentration 

L 

electrode  thickness 

Li+ 

lithium  ion 

r,  0 , 0 

spherical  coordinates 

R 

radius  of  spherical  particle 

S 

vacant  site  in  host  materials 

u 

radial  displacement 

X 

dimensionless  radial  position 

Greek  symbols 

Sr 

radial  strain 

£o 

tangential  strain 

0i 

volume  fraction  of  active  material  in  porous  elec¬ 
trode 

Or 

radial  stress 

°0 

tangential  stress 

o 

mean  stress 

CS  r,  max 

maximum  radial  stress 

^.max 

maximum  tangential  stress 

cr  shear,  max  maximum  shear  stress 

V 

Poisson’s  ratio 

Q 

partial  molar  volume  of  the  solute 

% 

dimensionless  stress  in  the  radial  direction  for 
potentiostatic  operation 

% 

dimensionless  stress  in  the  tangential  direction  for 
potentiostatic  operation 

4 

dimensionless  stress  in  the  radial  direction  for  gal¬ 
vanostatic  operation 

% 

dimensionless  stress  in  the  tangential  direction  for 
galvanostatic  operation 

X 

dimensionless  time 

rr 

dimensionless  total  elastic  energy  in  the  sphere  for 
potentiostatic  operation 

IT 

dimensionless  total  elastic  energy  in  the  sphere  for 
galvanostatic  operation 

and  surface  modulus  were  shown  to  play  a  significant  role  in  the 
determination  of  the  stress  amplitude  and  distribution  for  small 
particles.  For  the  purposes  of  this  work,  we  shall  not  include  the 
surface  phenomena,  allowing  us  to  streamline  the  exposition  and 
focus  on  the  distinguishing  differences  between  potentiostatic  and 
galvanostatic  operation. 


2.  Analysis  and  results 

2.1.  Mechanics  of  composition  induced  stresses 


We  consider  stress  caused  by  diffusion  within  a  spherical  parti¬ 
cle  of  radius  R.  The  bulk  of  the  spherical  particle  is  assumed  to  be 
an  isotropic  linear  elastic  solid.  Using  the  analogy  between  ther¬ 
mal  and  DIS  [1-7,17-20,22,23,25],  the  stress-stain  relationships, 
expressed  in  the  spherical  coordinate  system,  for  the  radial  and 
tangential  components,  are 

Sr  =  |( CTr  -  2 VOq)+  ^QC,  £0  =  g[(l  -  v)Oq  -  VOr]  +  ^QC, 

(1) 


where  E  is  Young’s  modulus,  v  is  Poisson’s  ratio,  C  is  molar  con¬ 
centration,  and  Q  is  the  partial  molar  volume  of  the  solute.  We 
further  assume  that  the  elastic  properties  are  independent  of  the 
concentration  C. 

Because  of  spherical  symmetry,  the  radial  and  tangential  strains, 
in  the  infinitesimal  formulation  of  deformation,  are  given  by 
£r  =  duldr  and  se  =  ulr,  where  u  is  the  radial  displacement.  Since 
atomic  diffusion  in  solids  is  a  much  slower  process  than  elastic 
deformation,  mechanical  equilibrium  is  established  much  faster 
than  that  of  diffusion.  Mechanical  equilibrium  is,  therefore,  treated 
as  a  static  equilibrium  problem.  In  the  absence  of  any  body-force, 
the  equation  for  static  mechanical  equilibrium  in  the  bulk  of  a 
sphere  is  given  by  Ref.  [25]: 


d(fr  ,  nGr-CTQ 
dr  r 


=  0. 


(2) 


The  solutions  for  the  normal  and  tangential  stresses  that  satisfy  the 
boundary  condition  or(R)  =  0  and  remain  finite  at  r = 0  are  given  by: 

2  FQ 

Vr(r)=^—V)lCav(R)-Cav(r)l 

FQ 

cre(r)  =  ^^[2 CaV(R)  +  Cay(r)  -  3C(r)],  (3) 


where  Cav(r)  =  (3/r3)  J^r2C(r)dr  is  the  average  concentration  in 
the  spherical  volume  of  radius  r  within  the  particle  of  radius  R.  Since 
limCav(r)  =  C(0),  Eq.  (3)  shows  that  the  stress-state  at  the  center  of 

r^O 

the  sphere  is  purely  hydrostatic,  i.e.,  lim<jr(r)  =  lim cr^r),  which  is 

r^O  r^O 

true  for  any  physically  admissible  concentration  profile. 

Recognizing  that  o0(r)  =  o^r),  we  can  write  the  “mean”  stress  as 


a(r)  =  gr(r)  +  2gg(r) 


2EQ 
9(1  -^) 


[Ca„(R)-C(r)]. 


(4) 


Because  of  the  spherical  symmetry,  one  principal  shear  stress  is 
zero  and  the  other  two  are  both  equal  to  (or  -  cre)l 2.  Eq.  (3)  shows 
that  the  principal  shear  stress  is 


or(r)-ae(r) 

2 


FQ 

w^)[c{r)-CM 


(5) 


Hence,  the  stresses  at  any  given  location  and  time  can  be  obtained 
once  composition  profile  is  known. 

From  the  stresses,  we  can  calculate  the  strain  energy  per  unit 
volume  or  strain  energy  density  e(r)  accumulated  as  a  result  of  the 
elastic  deformation  for  the  isotropically  deformed  sphere  [25]: 


a?{.r)  +  2aj{r)-2vae(r)[2ar(r)  +  ae{.r)] 
e(r)  = - - 2! - '  (6) 

The  total  elastic  energy  stored  in  the  sphere,  ET,  which  provides 
the  driving  force  for  fracture,  can  then  be  obtained  by  integrating 
the  strain  energy  density  over  the  entire  volume  of  the  spherical 
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particle: 

Ej  =  An  f  e{r)r2  dr.  (7) 

Jo 

2.2.  Diffusion-induced  stresses  under  constant  voltage  and 
current  conditions 

We  consider  reactions  of  the  type 

Li+  +  e”  +  S  =  [Li-S]  and  H+  +  e"  +  S  =  [H-S], 

reflecting  the  insertion  and  extraction  of  lithium  ions  (as  in  lithium 
ion  cells)  or  protons  (as  in  nickel  metal  hydride  cells),  respectively, 
two  cases  of  immediate  interest  for  battery  systems.  A  vacant  site 
with  the  host  materials  is  represented  by  S  and  e-  refers  to  the  elec¬ 
tron.  The  number  of  electrons  taking  place  in  the  electrochemical 
reaction  and  the  stoichiometric  coefficient  of  the  reacting  species 
are  taken  to  be  unity,  consistent  with  the  stated  elementary  reac¬ 
tions.  When  the  current  is  controlled  and  the  electrode  resistance 
is  governed  by  diffusion  limitations  within  the  insertion  particles, 
the  cell  current  IceU  is  related  to  the  current  over  a  spherical  particle 
by  I  =  Icenl{asL),  where  the  specific  surface  as  (commonly  expressed 
in  units  of  cm2  cm-3)  is  given  by  36>i/R  [26],  with  6\  representing 
the  volume  fraction  of  the  active  material  in  the  porous  electrode 
and  L  being  the  porous  electrode  thickness.  We  shall  assume  the 
current  density  is  uniform  over  the  particle  surface  and  we  shall 
not  account  for  changes  in  surface  area  of  the  particle  during  inser¬ 
tion  and  extraction  in  terms  of  impacting  the  current  density  [18], 
consistent  with  the  infinitesimal  formulation  of  deformation.  When 
the  electrode  voltage  is  controlled  relative  to  a  reference  electrode 
and,  again,  particle  resistance  dominates,  the  electrode  voltage  sets 
the  solute  surface  concentration  if  the  magnitude  of  the  dimension¬ 
less  group  i0  refRl[FD{CR  -  C0)]  is  large,  as  this  group  characterizes 
solute  diffusion  resistance  within  the  host  particle  to  that  of  the 
electrochemical  reaction  resistance  over  the  particle  surface  [26]. 
In  this  relation,  C0  is  the  initial  (uniform)  solute  concentration,  CR  is 
the  surface  concentration  during  the  constant  potential  operation, 
D  refers  to  the  diffusion  coefficient  of  the  solute,  and  F  is  Faraday’s 
constant.  The  exchange  current,  which  houses  the  interfacial  rate 
constants,  is  rendered  as  i0rej,  where  the  subscript  “ref’  indicates 
that  it  is  based  on  reference  concentrations.  General  texts  can  be 


ficient.  The  relevant  diffusion  equation  in  the  spherical  coordinate 
system  is  given  by  Ref.  [32]: 


dc  _  d  d_  f  2ac\ 

3 1  -  r2  dr  y  dr  J  ’ 


(8) 


The  initial  and  boundary  conditions  corresponding  to  potentiostatic 
control  are  given  by: 

C(r,  0)  =  C0,  forO  <r  <R 

C(R,  t)  =  CR,  for  t  >  0  .  (9) 

C(0,  t)  =  finite,  fort>0 


The  initial  and  boundary  conditions  corresponding  to  galvanostatic 
control  are  given  by: 

C(r,  0)  =  C0,  forO  <  r  <R 


D 


D 


d C(r,  t) 


dr 

d C(r,  t) 


dr 


=  1-,  fort  >  0 
F 

=  0,  fort  >  0 


(10) 


r=R 


2.2.1  Potentiostatic  (constant  voltage  and  surface  concentration ) 
operation 

The  analytic  solution  of  the  diffusion  equation  (Eq.  (8))  with  ini¬ 
tial  and  boundary  conditions  for  the  potentiostatic  control  (Eq.  (9)) 
is  [32]: 


C(r,  t)  -  C0 
Cr-Cq 


(—1  )n 

=  1  +  2 S'— —  sin(n7rx)e-n27r2\ 
2-^  nnx 

n= 1 


(11) 


where  x  =  r/R  and  r  =  Dt/R2 . 

The  average  concentrations  Cav(r)  and  Cav(R)  are  then  given  by 

OO  n 

Ca'J9  C°  =  i  +  [sin(mrx)  -  (n^x)cos(n^x)]e_n2jr2r 

Cr  -  Co  x6  ^  nini 

n=l 


and 


Cav(R)  ~  C0 
Cr-C0 


OO 

=  l-6y2^e-nVt. 


(12) 


Eq.  (12)  can  be  substituted  into  Eq.  (3)  to  obtain  the  solutions  for 
stresses  under  potentiostatic  control: 


Or 

[Etf(C*-C0)/(3(l-v))] 

_ <*e _ 

[E£2(Cr  -  C0)/(3(l  -  v))] 


4^^e-n27r2r 


1  f  sin(n7rx)  -  nnxcos{nnx)  \ 


(mry 

2 

( rwy 


{nnx) 


=  r) 


-2Ve-"2l2t  (  -2—  +  tlT  sin(njrx)  -  (-1)"  (  ^fmx)-mixcos(mtx)\ 
^  \  (nW  nnx  \  (njrx)3  J 


$(x,  r) 


(13) 


consulted  for  a  description  of  the  interfacial  charge-transfer  the¬ 
ory  [27,28].  To  summarize,  when  solute-diffusion  control  prevails 
and  the  electrode  is  subjected  to  galvanostatic  (controlled  current) 
operation,  the  current  density  over  the  particle  surface  /  =  FD(dC/dr) 
is  fixed.  Similarly,  when  solute-diffusion  control  prevails  and  the 
electrode  is  subjected  to  potentiostatic  (controlled  potential)  oper¬ 
ation  and  the  dimensionless  group  i0  re/R/[FD(Q  -  C0)]  is  large  (e.g., 
for  facile  interfacial  reactions,  electrodes  constructed  with  large 
particles,  a  low  solid-phase  lithium  diffusion  coefficient,  or  for 
any  combinations  of  these  situations),  the  surface  concentration 
is  a  constant  and  can  be  determined  from  an  equilibrium  poten¬ 
tial  relation.  A  specific  treatment  of  equilibrium  relations  based 
on  molecular  thermodynamics  for  intercalation  electrodes  can  be 
found  in  Refs.  [29-31  ]. 

To  illustrate  the  differences  in  stress  evolution  under  galvano¬ 
static  and  potentiostatic  control,  we  consider  a  simple  problem  of 
diffusion  within  a  sphere  of  radius  R  with  a  constant  diffusion  coef- 


where  f^(x,  r)  and  f^(x,  r)  are  the  respective  dimensionless 
stresses  in  the  radial  and  tangential  directions  for  potentiostatic 
operation. 

The  mean  and  shear  stresses  are  then  given  by 

(oy  +  2ae)/3 
[£^(Cr-C0)/(3(1-v))] 


 4’yy-n27T2T 

n=l 


(-1)" 

3nnx 


sin  (nnx)  j 


(14) 


and 


_ °r  —  G0 _ 

[2EQ(Cr  -  C0)/(3(l  -  v))] 


t  -i  yi 

El"1)  c-n27T2r 

nnx 

n= 1 


^sin(n7rx)  -  3 


sin(nnx)  -  nnxcos(nnx)\ 
{nnx)2  J 


(15) 
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Fig.  1.  Concentration  profile  (a)  and  the  corresponding  radial  (b),  tangential  (c),  and  shear  stress  (d)  for  potentiostatic  (constant  solute  surface  concentration)  control.  The 
solution  concentration  increases  with  time  (a),  and  the  magnitudes  of  the  tangential  and  shear  stresses  (c  and  d,  respectively)  at  the  particle  surface  decrease  with  time  from 
1  and  1/2,  respectively,  (cf.  Eq.  (3))  throughout  the  insertion  process. 


respectively.  Fig.  1  shows  the  concentration,  radial,  tangential,  and 
shear  stresses,  normalized  by  E£2(Cr  -  Co)/(3(l  -  v)),  as  a  function 
of  position  and  time.  Two  cases  of  special  interest  are  ( 1 )  “insertion” 
with  initial  condition  C0  =  0  and  a  finite  Cr  and  (2)  “extraction”  with 
a  finite  C0  and  CR  =  0.  For  insertion,  the  radial  stress  is  tensile  in  the 
sphere.  At  any  given  time,  the  tensile  radial  stress  is  the  highest 
at  the  center  and  decreases  monotonically  to  zero  at  the  surface 
(Fig.  lb).  It  can  be  shown  [25]  that  when  the  dimensionless  time 
r  =  Dt/R2  =  0.0574,  which  is  just  after  the  solute  reaches  the  center 
of  the  sphere,  the  radial  stress  at  the  center  attains  the  maximum: 


<77,  max 


E^(Cr-C0) 
3(1 -v) 


(16) 


The  radial  stress  at  the  center  then  deceases  with  increasing  time. 
The  tangential  stress  is  compressive  at  the  surface  and  tensile  at 
the  center  of  the  sphere  (Fig.  lc).  The  tensile  tangential  stress  at 
the  center  appears  before  the  solute  reaches  there.  At  the  center, 
the  tangential  and  radial  stresses  always  have  the  same  magnitude 
so  that  the  stress  at  the  center  is  purely  hydrostatic  in  tension.  The 
maximum  tangential  stress,  which  is  compressive,  occurs  at  the 
surface  at  time  zero: 


^max  — 


E£2{Cr  -  C0) 
3(1  -v) 


(17) 


Analogous  conclusions  hold  for  extraction  if  “tensile”  is  replaced 
by  “compressive”.  The  shear  stress  at  the  center  of  the  sphere  is 


zero.  For  all  times,  the  shear  stress  increases  towards  the  surface. 
The  maximum  shear  stress,  with  its  magnitude  half  of  that  of  the 
maximum  tangential  stress,  occurs  at  surface  at  time  zero  (Fig.  Id): 


EQ{Cr-Cq) 

O shear,  max-  6(1  _  y) 


(18) 


Knowing  the  stresses,  we  can  obtain  the  total  elastic  energy 
stored  in  the  sphere  ET  as  a  function  of  dimensionless  time  r  using 
Eqs.  (6),  (7)  and  (13): 


_ EtM _ 

27tR3E[£2(Cr  -  C0)/(3(l  -  v))]2 


-S' 


r)+2^"(x,  t)-2v$(x,  r)(2frv(x,  r)+^(x,  r)) 


fV2 


x2  dx 


=  nv(v,  t) 


(19) 


where  jqV  is  defined  as  the  dimensionless  total  elastic  energy  in 
the  sphere  for  potentiostatic  operation.  Eq.  (19)  is  shown  in  Fig.  2 
for  representative  values  of  Poisson’s  ratio  (-1  <  v  <  0.5)  for  poten- 
tiostatically  controlled  charging.  In  a  uniaxial  tensile-strain  test, 
the  Poisson  ratio  v  is  the  ratio  of  lateral  to  longitudinal  strain.  It  is 
dimensionless  and  typically  ranges  from  0.2  to  0.5.  For  incompress¬ 
ible  materials,  v  =  0.5,  and  the  material’s  volume  remains  constant 
under  deformation.  When  stretching  a  material  under  uniaxial  ten¬ 
sion  causes  no  lateral  contraction,  v  =  0.  For  any  stable  material, 
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0.0  0.1  0.2  0.3  0.4 


_ <*e _ 

[(£«/(3(l  -v)))(«/(fD))] 


=  i(l-2x2)  +  2]T 

n= 1 


exp(-A„r) 
X„  sin(A.„) 


(  sin(A,nx) 
\  Anx 


sin(A,nx)  -  (Anx)cos(A,nx) 

ip 


=  *)■ 


(22) 


where  £f(x,  r)  and  ^(x,  r)  are  the  respective  dimensionless  stress 
in  the  radial  and  tangential  directions  for  galvanostatic  operation. 
The  mean  and  shear  stresses  are  then  given  by 

{<*r  +  2oe)/3 

[(F£?/(3(l  -  v))){IR/{FD))] 


=  l(l-2x2)  + 


OO  / 

4  ^  /  sin(A,nx) 
3x“f  \  sin(A„) 


(23) 


Fig.  2.  Total  strain  energy  vs.  time  for  potentiostatic  charging.  The  strain  energy 
increases  monotonically  with  decreasing  Poisson  ratios. 


-1  <  v<0.5.  Poisson’s  ratio  affects  the  total  stored  energy  signifi¬ 
cantly.  The  total  strain  energy  first  increases,  reaches  a  maximum, 
and  then  decreases  because  of  the  transient  nature  of  the  stresses 
(Fig.lb-d). 


2.2.2.  Galvanostatic  (constant  current  and  surface  flux)  operation 
The  analytic  solution  of  the  diffusion  equation  (Eq.  (8))  with 
the  initial  and  boundary  conditions  (Eq.  (10))  for  the  galvanostatic 
condition  is  also  well  known  [32 J1 : 


C(r,  t) 


„  IR 
~Co  +  FD 


1  2  3 

3r  +  -x  -  — - > 

2  10  X2L^ 

n= 1 


sin(A,nx) 

A2sin(A.„) 


exp(-X2r) 


(20) 


where  Xn  (n  =  1,  2, 3, . . .)  are  the  positive  roots  of  tan(An)  =  Xn. 
Using  Eq.  (20),  we  can  express  the  average  concentrations  as 


Qy(r,  t) 


C0  +  3 


T+io(x2 


1) 


OO  / 

2  (  sin(Anx)  -  (Anx)cos(Anx) 

V  sin(A,„) 


exp(-A.2r) 


Cav(R,  t ) 


c»+3(S)'-c»+4'- 


(21) 


Eq.  (21 )  can  be  substituted  into  Eq.  (3)  to  obtain  the  solutions  for 
stresses: 


Or 

KEG/m  -  v))X«/(fD))] 


OO  / 

=  hi  y21  I  4  I  sin(^nX)  -  (A.„x)cos(A.nx) 

“5  x3trv  ^sin(A„) 


exp(-X2r) 


=  tf(x,  r), 


1  The  concentration  C(r,  t )  cannot  exceed  the  concentration  of  available  sites  for 
the  solute  in  the  case  of  insertion  electrodes,  and  this  will  place  a  limit  on  the 
admissible  value  of  the  product  It. 


and 


or  -  oe 

[(2£i2/(3(l  -  v)))(JR/(FD))] 


x 2  |  3  Tp  (  sin(Xnx)  -  (A.nx)cos(Xnx) 

10  +  ^2-.  (  Xpnp) 


exp(-A2r) 


OO  / 

1  ^  [  sin(A.nx) 
\A.2sin(A.n) 

n=l  x 


exp(-X2r) 


(24) 


respectively. 

Fig.  3  shows  the  concentration,  radial,  and  tangential  stresses, 
normalized  by  E£2IRI{ 3(  1  -  v)FD),  as  a  function  of  position  and  time 
and  using  the  sign  convention  defined  in  Eq.  (10)  for  the  current; 
i.e.,  insertion  and  extraction  correspond  to  a  positive  and  negative 
current,  respectively.  From  a  qualitative  perspective,  many  of  the 
observations  described  for  the  potentiostatic  case  (constant  surface 
concentration)  are  seen  for  the  galvanostatic  (constant  surface  flux) 
case  with  some  significant  exceptions. 

During  insertion,  the  radial  stress  is  tensile  in  the  sphere.  At  any 
given  time,  the  tensile  radial  stress  is  the  highest  at  the  center  and 
decreases  monotonically  to  zero  at  the  surface.  At  any  location,  the 
tensile  radial  stress  increases  with  time  and  tends  to  a  “steady- 
state”  (cf.  Eq.  (22)): 

Or,t^  oo  =  |(1  -X2) 


E£2 


(M) 


(25) 


The  maximum  tensile  stress  at  the  center  of  the  sphere  is 
o r.max  =  E*f2/F/(15(1  -  v)FD).  This  is  in  contrast  to  the  potentiostatic 
case  where  crr, max  exhibits  a  transient  behavior.  The  tangential 
stress  is  compressive  at  the  surface  and  tensile  at  the  center  of  the 
sphere.  The  tensile  tangential  stress  at  the  center  appears  before  the 
solute  reaches  there.  At  the  center,  the  tangential  and  radial  stresses 
always  have  the  same  magnitude  so  that  the  stress  at  the  center  is 
purely  hydrostatic  in  tension.  At  any  location,  the  tangential  stress 
increases  in  magnitude  with  time  and  tends  to  a  steady-state, 


^max  — 


1(1 -2x2) 


E£2  / IR\ 

3(1 -v)  \FDJ 


(26) 


The  maximum  tangential  stress  is  tensile  at  the  center  and  compres¬ 
sive  at  the  surface  with  the  same  magnitude  as  that  of  the  maximum 
radial  stress  at  the  center,  |F£?/F/(15(1  -  v)FD)|.  The  crossover  from 
tensile  to  compressive  tangential  stress  occurs  at  x  =  1  /V2  after 
the  stress  reaches  the  steady-state.  Similar  conclusions  (i.e.,  replac¬ 
ing  tension  by  compression)  hold  for  extraction.  The  shear  stress 
at  the  center  of  the  sphere  is  zero.  As  was  observed  for  the  poten¬ 
tiostatic  condition,  the  shear  stress  increases  towards  the  surface 
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Fig.  3.  Concentration  profile  (a)  and  the  corresponding  radial  (b),  tangential  (c),  and  shear  stress  (d)  for  galvanostatic  operation.  The  solution  concentration  increases  with 
time  (a),  and  the  magnitudes  of  the  tangential  and  shear  stresses  (c  and  d,  respectively)  at  the  particle  surface  increase  with  time  throughout  the  insertion  process.  By 
comparing  the  analogous  plots  of  Figs.  1  and  3,  one  can  discern  the  very  different  character  of  the  stress  evolution  for  potentiostatic  and  galvanostatic  control. 


for  all  times.  In  addition,  for  all  positions  within  the  sphere,  the 
shear  stress  increases  with  time.  The  maximum  shear  stress,  with 
its  magnitude  half  of  that  of  the  maximum  tangential  stress,  occurs 
at  surface  after  reaching  the  steady-state: 

1 

&  shear,  max  — 


EQ 

CM 


(m) 


(27) 


The  total  dimensionless  elastic  energy  stored  in  the  sphere,  77',  for 
galvanostatic  operation  can  be  expressed  as 


Et(  t) 


2jiR3E[f2IR/(3(\  -v)FD)] 

r>l 


f  &2(x,  T)+2^2 (x,  r)-2vf'(x,  t)(24(x,  r)  +  ^(x,  r)) 

Jo  L 


=  77'(v,  x) 


x 1  dx 


(28) 


Eq.  (28)  is  shown  in  Fig.  4  for  representative  values  of  Pois¬ 
son’s  ratio  for  galvanostatically  controlled  charging.  Similar  to  the 
potentiostatic  case,  Poisson’s  ratio  affects  the  total  stored  energy 
significantly.  Unlike  the  potentiostatic  case,  however,  the  total 
strain  energy  increases  monotonically  to  a  steady-state  consistent 
with  the  continuous  increase  in  the  stresses  with  time  (Fig.  3b-d). 


0.0  0.1  0.2  0.3  0.4 

T 


Fig.  4.  Total  strain  energy  vs.  time  for  galvanostatic  charging.  While  the  strain  energy 
increases  monotonically  with  decreasing  Poisson  ratios,  as  is  observed  in  the  poten¬ 
tiostatic  case  (cf.  Fig.  2),  the  evolution  of  the  strain  energy  during  galvanostatic 
control  is  quite  different  from  that  of  potentiostatic  control. 
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2.3.  Implications  for  cell  operation 


There  have  been  several  recent  studies  optimized  charging  con¬ 
ditions  [33-35].  Many  of  these  studies  focus  on  shortening  the  time 
to  achieve  a  high  state  of  charge  (e.g.,  Ref.  [33]).  There  is,  however, 
little  discussion  on  optimizing  charging  conditions  to  minimize 
mechanical  stresses  within  the  electrodes.  The  above  analysis  pro¬ 
vides  a  guide  to  deciding  between  constant  current  and  constant 
potential  charging  methods.  Specifically,  we  can  compare  the  max¬ 
imum  stresses  encountered  during  either  purely  potentiostatic  or 
galvanostatic  conditions,  subject  to  the  assumptions  used  in  the 
formulation  the  equations  presented  previously.  The  ratio  of  max¬ 
imum  radial  stresses  at  the  center  is  given  by 


potentiostatic,  max 


galvanostatic,  max 


nCR-C0 
IR/(FD ) 


(29) 


Thus,  the  maximum  radial  stress  under  potentiostatic  condition 
is  greater  than  that  under  galvanostatic  condition  if  the  magni¬ 
tude  of  current  I  is  less  than  2FD(CR  -  C0)/R  and  vice  versa.  Eq.  (29) 
corresponds  to  the  maximum  radial  stress  shown  in  Fig.  lb  (near 
r  =  0.0574,  cf.  Eq.  (17))  divided  by  the  maximum  radial  stress  of 
Fig.  3b  (obtained  from  the  steady-state  solution,  E£2IR/(  15(1  -  v)FD) 
per  Eq.  (25)).  The  constant  2  can  be  replaced  by  1.93  if  3-digit  accu¬ 
racy  is  desired. 

The  ratios  of  maximum  tangential  stresses  and  shear  stresses  at 
the  surface  are  the  same  and  are  given  by 


potentiostatic  potentiostatic 

°_0 _  =  °shear _  =  5  CR  ~  C0 

galvanostatic  galvanostatic  IR/(FD ) 

° shear 


(30) 


Thus,  the  maximum  tangential  or  shear  stress  under  potentiostatic 
condition  is  greater  than  that  under  galvanostatic  condition  if  the 
magnitude  of  the  current  I  is  less  than  5  FD(  CR  -  C0  )/R  and  vice  versa. 
Eq.  (30)  corresponds  to  the  maximum  tangential  and  shear  stresses 
for  a  constant  surface  concentration  (1  and  1/2,  respectively,  as  can 
be  seen  by  substituting  C0  for  Cav{R)  and  CR  for  C(R)  in  Eq.  (3)  for 
t= 0),  respectively,  divided  by  the  maximum  tangential  and  shear 
stresses  of  Fig.  3c  and  d,  corresponding  to  EZ2IRK  15(  1  -  v)FD)  for  the 
tangential  stress  and  half  this  value  for  the  shear  stress,  as  described 
below  Eq.  (26). 

In  Fig.  5,  we  plot  Eqs.  (29)  and  (30)  and  show  that 
(Cfl  -  C0)FD/(/R)  uniquely  determines  the  maximum  stress  ratios 
for  the  radial,  tangential,  and  shear  stresses.  The  horizontal  line 
at  (jpotentiostatic ^galvanostatic  =  ^  delineates  the  parameter  range  of 
(CR  -  C0)FD/(/R )  such  that  ^potentiostatic ^galvanostatic  is  greater  than 


or  less  than  1.  Similarly,  the  ratio  of  maximum  total  strain  energy 
can  be  obtained  using  Eqs.  (19)  and  (28).  Hence,  the  maxima  in 
the  curves  depicted  in  Fig.  2  relative  to  those  plotted  of  Fig.  4  (the 
long-time  solution  in  the  case  of  Fig.  4)  yields  the  points  depicted 
in  Fig.  6  for  [(IR/FD)/(Cr  -  C0)]2(Ff tentiostatic^ 

This  ratio  of  the  maxima  in  total  strain  is  seen  to  range  from  4  to 
5.5  for  admissible  values  of  the  Poisson  ratio. 

These  results  show  that  {CR-  Co)FD/(JR)  uniquely  determines 
whether  potentiostatic  or  galvanostatic  charging  would  gener¬ 
ate  greater  radial,  tangential,  and  shear  stress.  It  also  determines 
which  method  generates  greater  total  strain  energy  that  can  drive 
fracture  in  the  spherical  particle.  This  dimensionless  parameter 
(Cr  -  C0)FD/(/R),  representing  the  combined  effects  of  materials 
properties  (e.g.,  D),  size  (e.g.,  R),  and  charging  conditions  (e.g.,  CR ,  C0 
and  /),  can  therefore  be  used  to  guide  the  selection  of  charging  con¬ 
ditions  to  minimize  the  maximum  stresses  and  total  strain  energy 
experienced  by  an  insertion  electrode.  We  consider,  for  example, 
an  electrode  made  of  a  brittle  ceramic  material  that  fractures  due 
to  tensile  stresses.  Eq.  (29)  can  then  be  used  to  select  a  charging 
condition,  either  purely  potentiostatic  or  galvanostatic,  which  gives 
rise  to  a  lower  tensile  stress.  In  contrast,  since  plastic  deformation 


Fig.  5.  A  map  showing  the  ratio  of  maximum  stresses  under  potentiostatic  and  gal¬ 
vanostatic  control  vs.  a  dimensionless  parameter,  (Cr  -  Co)FD/(IR),  which  represents 
the  combined  effects  of  materials  properties  (e.g.,  D),  size  (e.g.,  R ),  and  charging  con¬ 
ditions  (e.g.,  CR,  Co,  and  /).  Hence  this  map  allows  one  to  comprehend  how  changes 
material  or  process  parameters  can  be  employed  to  stay  below  specified  stress  val¬ 
ues.  The  constant  2  appearing  in  the  label  for  the  radial  stress  ratio  can  be  replaced 
by  1.93  if  3-digit  accuracy  is  desired. 


Fig.  6.  Ratio  of  maximum  total  strain  energy  as  determined  by  Eqs.  (19)  and  (28). 
The  dashed  line  is  provided  as  a  guide  and  does  not  represent  calculations. 


is  usually  the  dominant  failure  mechanism  for  metallic  alloy  elec¬ 
trodes,  Eq.  (30)  can  be  used  to  select  a  charging  condition  resulting 
in  a  lower  shear  stress.  Thus,  depending  on  the  failure  mode  of 
the  electrode  material,  e.g.,  tensile  fracture  or  shear  induced  plastic 
deformation,  an  appropriate  charging  method  may  be  selected,  and 
the  best  regime  may  consist  of  a  combination  of  potentiostatic  and 
galvanostatic  conditions. 


3.  Summary 

We  have  developed  analytic  expressions  for  stress  evolution  in 
a  spherically  shaped  electrode  element  under  either  galvanostatic 
or  potentiostatic  conditions.  We  show  that  stresses  evolve  quite 
differently  for  the  two  conditions  when  solute  diffusion  resistance 
dominates  within  host  particles. 

•  For  potentiostatic  control,  the  magnitude  of  radial  and  tangen¬ 
tial  stresses  at  any  location  increases  initially  and  then  decreases 
with  time.  The  total  strain  energy,  which  drives  particle  fracture, 
experiences  a  similar  trend. 
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•  For  galvanostatic  control,  the  magnitude  of  radial  and  tangential 
stresses  at  any  location  increases  and  then  reaches  a  steady-state 
value  that  is  independent  of  time.  The  total  strain  energy  also 
increases  monotonically  to  a  finite  value. 

The  maximum  stresses  and  total  strain  energy  attained  under 
potentiostatic  and  galvanostatic  conditions  are  compared  and  con¬ 
trasted.  This  analysis  suggests  a  possibility  of  developing  new 
battery  charging  strategies  that  minimize  stresses  and  thus  prolong 
battery  life. 
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